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We have measured the angular diameter of the transiting extrasolar planet 
host star HD 189733 using the CHARA O/IR interferometric array. Combining 
our new angular diameter of 0.377 ±0.024 mas with the Hipparcos parallax leads 
to a linear radius for the host star of 0.779 ± 0.052 i?s un and a radius for the 
planet of 1.19 ± 0.08 i?j up . Adopting the mass of the planet as derived by its 
discoverers, we derive a mean density of the planet of 0.91 ± 0.18 g cm -3 . This 
is the first determination of the diameter of an extrasolar planet through purely 
direct means. 

Subject headings: infrared: stars, stars: fundamental parameters, techniques: 
interferometric, stars: individual: HD 189733 



Introduction 



A handful of extrasolar planets transit their host stars, causing a reduction in stellar flux 
as the planet blocks part of the star's disk. The planet orbiting HD 189733 is one of the four- 
teen known transiting planets. Usi ng radial velocity an d photometric measurements made at 
the Haute-Provence Observatory, iBouchy et al.l (120051 ) discovered a hot Jupiter-like planet 
with an orbital period of 2.219 days and estimated the star's radius to be 0.76 ± 0.01 R Q . 
This value, along with a planet-to-s tar radius ratio of 0. 172 ± 0.003, led to a planetary radius 
of 1.26 ± 0.03 -Rjup- More recently, iBakos et ajj (l2006bl ) refined the orbital parameters using 
BVRI multi-band photometry and found the planet's radius to be 1.154 ± 0.032 i?j up . 

We observed HD 189733 using Georgia State University's Center for High Angular 
Resolution Astronomy (CHARA) Array in order to directly determine the host star's radius 
and thereby calculate, in a strictly geometric manner, the radius and density for the planet. 

Planetary densities were previously estimated from photo metric observations of the 
transiting planets and range from 0.38 k cm" 3 for HD 209458b (ICharbonneau et al.ll2000i ) 
to 1.17 g cm -3 for HD 149026b (ISato et al.l 120051 ). These density calculations are highly 
dependent on estimated stellar diameters based on spectral energy distribution (SED) fits 
using published photometric values, which are fundamentally indirect in nature, relying 
upon a priori assumptions regarding the host stars' stellar atmospheres. For the four 'bright' 
(V < 12) transit host stars, these angular sizes are in the range of 0.05 to 0.40 milliarcseconds 
(mas). The longest baselines of the CHARA Array are capable of resolving the largest and 
brightest of these objects. 
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2. Interferometric Observations and Diameter Determination 
2.1. Observations and Data Reduction 

Spatially resolved observations of HD 189733 were obtained w ith the CHARA Array . 



a six-element interferometer located on Mount Wilson, California (IMcAlister et al.ll2005l ). 
The Array operates in two wavelength regimes: in visible wavelengths (470-800 nm) for 
tracking and tip/tilt correction; and in the near infrared K' (2.13 /xm) and H (1.67 /im) 
bands for fringe detection. Because of the small angular diameter for the star, only the 
if -band observations obtained at our longest baseline pair (telescopes El and SI) are used 
in our final diameter analysis. 

HD 189733 was observed on several nights during the summer of 2006 along with the 
calibrator star HD 190993, a B3 V star offset by 1.7°, selected on the basis of its small 
estimated angular diameter and its apparent lack of any close companion. The latter cri- 
terion was verified by a thorough literature search while a spectral energy distribution fit 
to HD 190993 led to an estimated angular diameter of 0.167 ± 0.035 mas with no residuals 
suggestive of a companion (see Figured]). This results in a predicted visibility (V) for the 
calibrator of Kai=0.961+°;°^ at our longest baseline of 330 m, resulting in a contribution 
of oy ~ 0.01 — 0.02 to the calibrated visibility errors seen in Table 1. The small angular 
size and high visibility of the calibrator means HD 190993 is essentially unresolved using the 
CHARA Array, and the uncertainty in visibility due to calibrator diameter error is small 
compared to the measurement error. Therefore uncertainties i n the calibrator diameter wi ll 



not affect the HD 189733 diameter measurement significantly (Ivan Belle fc van Belle! 120051 ) . 
Even HD 190993's considerable v sin i does not contribute error to our diameter fits due to 
its small angular size. 



We note that the M-dwarf companion to HD 189733 reported by iBakos et al.l (l2006al ) 
on the basis of common space motion at an angular separation of 11.2 arcsec is well outside 
the interferometric field of view, and its presence has no effect on our results. Although the 
effect on visibility would be small in the first lobe of the V(B) curve, we have confirmed that 
our observed epochs do not occur within the predicted times of planetary tr ansit or eclipse 



using the period and reference time of central transit of IBakos et al.l (j2006bl ) 



All our observations were obtained with the single-baseline, pupil-plane "CHARA Clas- 
sic" beam combiner, and we employed the standard practice of observing the target and 
calibrator sequentially to provide a series of time-bracketed observations from which the 
instrumental visibilities could be reduced to calibrated values for the target star. The 
observing practice and redu ction process employed here is identical to that described by 



ten Brummelaar et al.l (120051 ). The results of this process are summarized in Table [TJ 
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Single baseline Michelson stellar interferometers measure complex visibilities, usually 
recorded as amplitudes and phases, which are related to the intensity distribution of the tar- 
get through a Fourier transform; phase information is typic ally corrupted by the atmosphere, 
leaving the amplitude, referred to simply as the visibility (jColavita et al.l Il999f ) . 



2.2. Diameter Fit 



Diameter fits to the visibilities and baselines from Table [T] were performed using the 
uniform disk (UD) approximation given by: 



V 



2Ji(x) 



x 



where J\ is the first-order Bessel function and 

x = 7t5^udA _1 , 



(2) 



where B is the projected baseline at the star's position, #ud is the apparent UD angular 
diameter of the star, and A is the wavelength of the observation. T he limb-darkened (LP) 



relat ionship incorporating the linear limb darkening coefficient (IHanbury Brown et al. 
19741 ) is given by: 
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These fits resulted in 6 UD = 0.366 ± 0.024 mas a n d 8 Tr > = 0.377 ± 0.024 mas, the latter 
incorporating fix = 0.36 taken from I Claret et al.l (119951 ) after adopting log g = 4.5 and 
T e ff = 5000 K for HD 189733 (see Figure [2]). The reduced x 2 minimization in both cases 
yielded a value of 1.593, and the errors quoted are for an increase of the x 2 value of 1.0, that 
is, the 68% confidence interval. Dividing this x 2 by the number of degrees of freedom, which 
in our case is 8, yields 0.199, which is much less than 1.0 showing that the fit is quite good 
and that our error estimates for the visibility points are conservative. If we rescale these 
errors bars to force x 2 to be equal to the number of degrees of freedom, which assumes that 
there are no systematics in the measurements, they are approximately half the size as they 
are shown in Figure 2 and would also reduce our final error estimates by a factor of two. 
However, we will remain conservative and continue to use the error estimate based on the 
raw x 2 value. 
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2.3. Estimate of the Angular Size of HD 189733 



An a priori estimate of the angular size of HD 189733 is a parameter of considerable 
interest, because the size of HD 189733b is de termined only relati ve to the size of its parent 
star from the photometric transit timing data. IBakos et al.l (j2006bl ) consider no less than four 
separate methods in t heir investigation of the system: V — K color angular radiu s predi ction 
(Kcrvella et al. temperature radius, isochrone radii from iGirardi et al.l (|2002i) an d 



Baraffe et al. 



(119981 ). and the Johnson V - 2MASS Teff calibration of iMasana et al.l (120061 ). 



None of these approaches appears to have much merit, sin ce the only p rimary data 
we ha ve been able to find in the literature were Tyc ho Br, Vt (IBesselll l2000l ) . Stromgren 



ubvy (folsenlll993h . and 2MASS JHK photometry JCutri et all l2003h . No spectroscopy 



or measures of log g appear to be availabl e in the literature, n or d o direct measures o f 
Johnson photometry. The values of V used in lKervella et al.l (120041 ) and lMasana et al.l (120061 ) 
appear to have been extrapolated from Vt- Furthermore, sustained long-term observations of 
HD 1897 33 by the asteroseismology satellite have found the star to be photometrically 

variable (lRowdl2007l ). casting significant doubt on any radius derived from a photometric 
relationship. Based on th is informatio n , we co nsider the size errors for HD 189733 quoted in 
Bouchy et al.l (120051 ) and IBakos et al.l (j2006al ) derived from the methods cited above to be 
underestimates. 

To explore what is a more appropriate error for an inferred angular size, we executed a 
spectral energy distribution fit of the available spectrophotometry for HD 189733 cited above. 
Given the known variability of HD 189733, the quoted millimagnitude error estimates of the 
Tycho and Stromgren photometric data points were increased by a factor of ten. These 
phot ometric data p oints were fit to the solar-abundance K0 V and K2 V templates available 
from iPickled (119981 ). with the resulting fit values for reddening, bolometric flux, and angular 
diameter seen in Table [2j along with the appropriate x 2 P er degree of freedom (x 2 PDF) 
values. Unfortunately, a Kl V template is not available in that library, although an estimate 
of one from interpolating between the two bracketing spectral types was synthesized by us 
for testing this spectral type. These fits are seen in Figure [3J 



The appropriate model spectrum from lMunari et al.l (120051 ) for a 5000 K star was fit with 



a x 2 PDF of 2.80, but this model (nor any other available in the literature) unfortunately 



only covered the visible portion of the spectrum. The 115 to 2500 nm range of iPicklea (119981 ) 
was necessary to fully characterize the available photometry, and thus we constrained our 
analysis to this particular set of stellar templates. 



Our finding is that, even with this highly detailed analysis of the stellar spectral energy 
distribution, the most appropriate modelling of that SED reveals a predicted angular size of 
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only 9 = 0.363 ± 0.011 mas - a 3% error bar - which corresponds to a stellar linear radius of 
R = 0.752 ± 0.026 R^. 



Discussion 



Our new direct determination for the angular diameter of HD 189733 of 0ld — 0.377 ± 
0.024 mas can be com bined with the Hipparcos parallax for the star of ir = 51.9 ± 0.9 mas 
(jPerryman et al.lll997l ) to give a physical radi us for the star of /L t ar = 0.779 ± 0.052 R Q , 
which is about 3% larger than that adopted by iBakos et al.l (j2006bl ). 



By the nature of the light curve analysis, the relative increase for the radius of the 
host star will directly translate into the sam e relative in c rease in the radius of the planet 
HD 189733b. Thus, revising the radius of IBakos et al.l (j2006bl ) of 1.154 i?j up , our new 
estimate for this va lue is -R p i ane t = 1-19 ± 0.08 -Rj up . Furthermore, adopting the value of 
Bouchy et al.l (120051 ) for the mass of the planet of 1.15 Mj up , we derive a new estimate for 
the density of HP 1897 33b to be p = 0.91 ±0.18 g cm -3 . These values are in good agreement 
with lWinn et al.l (120071 ). who used transit photometry to constrain the stellar and planetary 
radii. The values of M p i anet , -R p i ane t, and p p i ane t are all consistent with the modest collection 
of these parameters presently available for transiting exoplanet systems and support the 
conclusion that HD 189733 is not among the few hot Jupiters that present extraordinarily 
large radii for their masses. 



4. Interferometric Non-Detection of Binarity of HD 189733 



Given the higher resolution of interferometric arrays, a possible close-separation tertiary 
companion may affect our measures of HD 189733's visibility and thereby complicate our 
interpretation. As such, it was pru dent for us to also o bserve HD 189733 with the Palo- 
mar Testbed Interferometer (PTI) (jColavita et al.l Il999f ) . an instrument with intermediate 
baselines on a variety of sky projections, suitable for exploration of possible unseen nearby 
luminous (stellar) comp anions. PTI has be en demonstrated to be sensitive to nearby com- 
panions with AK < 4. (Boden et al.l 19981) , which for a K2-3 V primary star rules out any 
M-dwarf companions (IBessell fc Brettlll988l ). 



PTI observed HD 189733 in the if-band on the nights of 2006 June 10-12, 2006 June 24, 
and 2006 July 8-10. Four of those nights used PTI's 85-m NW baseline configuration, two 
used the 110-m NS baseline, and one night was a 85-m SW baseline night. For all of these 
nights, HD 189733's normalized V data points were indistinguishable from unit visibility, 
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which corresponds to a completely unresolved point source, as would be expected for a single 
~ 0.37 mas star being observed by PTI at 2.2 ^m. 

5. Conclusion 

Our results for the radii of the host star and planet in the HD 189733 exoplanet system 
are formally in good agreement with existing measurements of these parameters as well as 
with the estimate for the density of the planet and have the additional and significant merit 
that they represent direct measurements of stellar and planetary diameters that do not rely 
upon inferences about stellar atmospheres. While the diameter measurements are currently 
at a 6% level of accuracy, we expect to improve this considerably as we implement fringe 
detection at shorter wavelengths at the CHARA Array. In the meantime, these results 
demonstrate a new role that long-baseline optical/infrared interferometry can play in the 
field of exoplanet astronomy. 
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at CDS, Strasbourg, France, and of NASA's Astrophysics Data System. This publication 
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Institute of Technology, funded by the National Aeronautics and Space Administration and 
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Table 1. Interferometric Measurements of HD 189733. 



MJD 


Baseline 


Visibility 


aVis 


(53886.0 +) 


(m) 






0.905 


330.5 


0.851 


0.071 


1.936 


327.9 


0.843 


0.056 


1.958 


324.9 


0.857 


0.054 


8.865 


330.5 


0.869 


0.034 


76.742 


326.5 


0.909 


0.069 


76.761 


323.8 


0.863 


0.049 


76.778 


321.3 


0.877 


0.045 


76.793 


319.0 


0.839 


0.045 


76.824 


315.5 


0.829 


0.061 
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Table 2. Spectral energy distribution fits for HD 189733 photometry to empirical spectral 

templates. 

Model Parameters 81 Fitted parameters 

Spectral type Effective Temperature x 2 PDF Reddening Bolometric Flux -Fbol Angular Diameter 
Teff (K) Ay (mag) (erg cm" 2 s" 1 fim^ 1 ) 9 (mas) 



KO V 5188 ± 50 

Kl V b 5040 ± 70 

K2 V 4887 ± 50 



2.62 0.308 ±0.031 3.207 
1.94 0.105 ±0.031 2.828 
2.39 0.00 ±0.030 2.613 



± 0.061 x 10" 8 0.365 ± 0.008 

± 0.049 x 10" 8 0.363 ± 0.011 

± 0.043 x 10" 8 0.371 ± 0.008 



a Models from lPicklesI l|l998h . 

b Kl V model interpolated from the K0 V and K2 V models. 
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Fit x PDF: 0.56 ( 10.7 / 19 DF) 
Est F bo , : 1.219e-06 +/- 8.69e-09 erg cm" 1 s 
Model T eff : 19055 +/- 2000 K 

Est Diam : 0.167 +/- 0.035 mas 
Est A v : 0.206 +/- 0.0133 mag 



Template SED from ukb3v.dat 




1 1.5 
Wavelength (fim) 



Fig. 1. — SED fit of our calibrator star, HD 190993. In the top panel, vertical bars represent 
the errors for the data points they overlay and the horizonal bars represent the bandpass of 
the data point. In the bottom panel, the fractional residuals (difference between data point 
and fit point, normalized by that data point) are shown for each data point. 
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Fig. 2. — The % 2 fit to HD 189733's visibilities. The solid line represents a theoretical 
visibility curve for a star with a limb-darkened diameter of 0.376 mas, the boxes are the 
measured visibilities, and the vertical lines are the measured errors. 
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HD 189733 Net SED Model 
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Fit % PDF: 1.94 (13.6 /7 DF) 
EstF bd : 2.828e-08 +/-4.?le-10 erg c 
Model T ff : 5040+-/- 70 K 
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HD 189733 Net SED Model 




Input Photometry i 1 — 

Model Photometry X 
Model Net (Source + Trans) SED 



Fit x PDF: 2.39 (16-7 11 DF) 
EstF bd :2.613e-08+/-4.32e-10erf 
Model T eff : 4887 +-/- 50 K. 

iam: 0.371 +/- 0.00818 mas 
: 0.000 +/- 0.0326 mag 



Template SED from ukic2v.dat 



1 1.5 
Wavelength (Jim) 



Fig. 3. — SED fits of HD 189733. The top panel shows the photometric fit for a KO V, the 
middle panel shows the Kl V fit, and the bottom panel shows the K2 V fit. All three fits 
allowed the reddening factor to vary. 



